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Abstract. The successful quasi-particle model is compared with recent lattice data 
of the coefficients in the Taylor series expansion of the excess pressure at finite 
temperature and baryon density. A chain of approximations, starting from QCD to 
arrive at the model expressions for the entropy density, is presented. 
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1. Introduction 



The equation of state of strongly interacting matter at finite temperature T and small 
O ■ chemical potential /i has become accessible fairly recently through first principle lattice 
QCD calculations [TJ El H|- Apart from solving QCD on the lattice, there exist 
analytical approaches such as resummed HTL scheme, $ functional approach etc. (cf. [3] 
for a recent survey). Ab initio approaches |E1IZ| are restricted in describing lattice data 
on temperatures T > 2.5T C , where T c is the transition temperature of deconfinement 
and chiral symmetry restoration. In contrast, phenomenological models with adjustable 
parameters [HI El EH cover the region T > T c . Here, we present new developments of 
our quasi-particle model (QPM) [H Ell- 
in section [21 the QPM is reviewed. In section El the model is confronted with 
recent lattice QCD data. Motivating our model, section [3] briefly illustrates a chain of 
approximations within a $ funtional approach. The results are summarized in sectional 

2. Quasi-particle model 

The model is based on the idea that the quark-gluon fluid can be expressed in terms of 
quasi-particles. The pressure of N q light (q), strange (s) quarks and gluons (g) reads 

p(T,fi)= Pi(T,p)-B{T,ti). (1) 

i = q,s,g 

Pi are thermodynamic standard expressions in which T and fi dependent self-energies 
Hi enter. B(T,p), together with the stationarity condition 8p/8mf = [H], ensures 
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thermodynamic self-consistency (cf. for details). Thus, the entropy density 
s = dp/dT = J2i = Q s g s i explicitly reads 

di r dkk 2 | (Ik 2 + mf) [fm + f _ (k)] _ ^ [fm _ f _ (k)] 1 (2) 



2n 2 T J ' [ uJi(k) 

with statistical distribution functions fi(k) = (exp[(wj(fc)q=/^)/T] + S'j)~ 1 , S g = S s = 1, 
S'g = —1, fi q = p,, fig = fig = and degeneracies d g = 6iV 9 , c? s = 6 and d g = 8. 

In the thermodynamically relevant region of momenta k ~ T, ^, the quasi-particle 
dispersion relations are approximated by the asymptotic mass shell expressions near 
the light cone uf(k) = k 2 + m 2 with m 2 = Ui(k;T,fi) + (a^T) 2 . The self-energies Ilj 
of the quasi-particle excitations are approximated by their 1-loop expressions at hard 
momenta jU] neglecting imaginary parts, and x{F represent the quark masses used on the 
lattice pQ. Replacing the running coupling in the self-energies by an effective coupling, 
G 2 (T, /i), non-perturbative effects are thought to be taken into account [T2J. Imposing 
thermodynamic consistency onto p, a flow equation for G 2 (T,/i) follows [9 j 

dG 2 dG 2 , 

a ^ + aT ^T =b (3) 
which can be solved as Cauchy problem by knowing G 2 on an arbitrary curve T(/a). 
One convenient choice is parametrizing G 2 {T{^i = 0)) appropriately (cf. such that 
p and s can be computed at non- vanishing \x from (|1|2|) . 

3. Expansion coefficients 

Apart from (jljl . the pressure can be decomposed into a Taylor series 

T 4 = ^4 + ^ = Z^ C2 "( T HtJ ^ 

n=0 

with Co(T) = p(T,fi = 0)/T 4 and vanishing c& for odd k. The expansion coefficients 
Cfc have been subject of recent lattice evaluations |T] by computing derivatives of the 
thermodynamic potential. They follow from ((T} as Ck(T) = d h p/d/i k \^_ T k ~ 4 /k\ in our 
model and read 

oat rco uj q /T 

ci = Sk dfcJb a -g- , ^ 2 , (5) 
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POO p^q/T ( T 1 

/ dk A; 2 - A e 2 ^ T - 4e^ T + 1 - — [e 2 ^/ T - l\ A 2 , (6) 

47r 2 T3y (e<V^ + l) 4 I u q 1 / 

where A 2 = (G 2 /tt 2 + 3x q ^jG 2 /(6ii 4 ) + [\x q T 2 /VQG 2 + \T 2 ]d 2 G 2 / '<V)| M=0 , w <? is taken 
at = and d 2 G 2 /dn 2 follows from differentiating ©. 

In the left panel of Figure HJ the QPM results of 02,4 calculated from (|5|(jjl are 
compared with lattice QCD results pQ for the two-flavour case, i.e. N g = 2. Using 
x q = 0.4 and x g = as in [I] and setting T c (/x = 0) = T = 170 MeV, G 2 (T(/x = 0)) 
is adjusted to describe c^iT). Since C4 in (jOJ only depends on G 2 and its derivatives at 
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/i = 0, no further assumptions enter into the evaluation once the parametrization is fixed. 
A fairly good agreement is found for c±(T). Note, in particular, that the pronounced 
peak structure of C4 at T solely originates from the term including d 2 G 2 /dfi 2 \ . 
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Figure 1. Left: Expansion coefficients C2(T) (squares) and ca(T) (circles). 
Data from J,. Full and dashed curves represent, respectively, corresponding QPM 
results ()5I6JI . Right: Scaled excess pressure Ap from £Q (symbols) and QPM results 
(full lines) from truncating Q at order (fi/T) 4 for different [i/T = 0.2, 0.4, 0.6, 0.8, 1.0 
(lower to upper curve). Dashed lines depict corresponding full QPM results. 



In the excess pressure Ap has been calculated as the truncation of the 
expansion Q including the order (/z/T) 4 . The right panel of Figure ^ exhibits the 
comparison of Ap calculated by employing only 02,4 in (j3J) (full lines) with the lattice 
data for different /i/To. An impressively good agreement with the data is observed for 
small values of [i in which case C4 is of less importance. Similarly, Ap can be evaluated 
as infinite series from the QPM by combining (0) and (@J) (dashed lines). These full 
results differ noticeably from the truncated results only for T « To. It should be noted 
that the model is successfully applied to describing the equation of state with strange 
quarks [T4*j . 

4. Contact with QCD 

Motivating the strong assumptions made in formulating the QPM expressions (ll|2|) . 
a chain of reasonable approximations starting from QCD would be of desire. The 
thermodynamic potential Q = —pV in ghost free gauge reads 

Q[D,S] =T^Tr [In IT 1 -UD] — Tr [in S^ 1 — ES] j + T$[D, S] (7) 

with dressed propagators D and S of bosons and fermions and corresponding exact 
self-energies IT and £ from Dyson's equations. The functional $ is given by the sum 
over all 2 particle irreducible skeleton diagrams and IT, E follow from truncating dressed 
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propagator lines in these diagrams |B] . The sum over Matsubara frequencies in the trace 
Tr is transformed into an appropriate contour integral in the complex energy plane. 
Computing the entropy density s = —d(Q/V)/dT, an ultra-violet finite expression 



is derived. After truncating $ at 2 loop order, s' = is found. Lost gauge invariance 
gets restored by employing hard thermal loop expressions for IT and X which show the 
correct limiting behaviour for k ~ T, fi. Performing the remaining trace tr over discrete 
indices in (JHJ), quantum numbers of quarks and gluons are recovered. Furthermore, 
the exponentially suppressed longitudinal gluon modes and the plasmino branch can be 
neglected. In addition, neglecting imaginary parts in II and X as well as Landau damping 
and approximating self-energies and dispersion relations suitably, expression (J2J) for the 
entropy density s is recovered. 

5. Conclusion 

The quasi-particle model has been reviewed and successfully compared with recent 
lattice data of the expansion coefficients 02,4 and the excess pressure at finite temperature 
and chemical potential. Briefly, a chain of approximations within the $ functional 
scheme starting from QCD has been summarized which leads to our model. 
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